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ABSTRACT 

We study the metallicity properties of galaxies in the zCOSMOS sample between 
0.2 < z < 0.8. At z < 0.46, where Ha and [Nil] are detected, we find the same depen- 
dence of metallicity on stellar mass and Star For mati on Rate (SFR), the Fundamental 
Metallicity Relation, found by Mannucci et al. (2010) in SDSS galaxies on a similar 
redshift range. We extend this relation to higher redshift, 0.49 < z < 0.8 where the 
R23 metallicity index can be measured in our data, finding no evidence for evolution, 
and a metallicity scatter around the relation of about 0.16 dex. This result confirms, 
with a much higher level of significance with respect to previous works, the absence 
of evolution of the FMR during the last half of cosmic history. 
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1 INTRODUCTION 

The Fundamental Metallicity Relation (FMR) is a tight 
relation between stellar mass , SFR , and metallicity, dis- 
covered by Mannucci et al. l|2010l . hereafter MIO) using 
local SDSS galaxies. The gas-phase metallicity of ~140.000 
galajdes with redshift between 0.07 and 0.3 is found to 
depend not only on stellar mass, i n the well-known mass- 
metallicity relation (Tremonti et al. l2004h . but also on SFR, 
with more active galaxies showing lower metallicities. After 
considering this relation, the residual metallicity scatter is 
limited to ~0.05 dex, which is compatible with the effect of 
the uncertainties on mass, SFR, and metallicity. 
The role of specific SFR (SSFR=SFR/Mass)) in r elatio n 
with metallicity was first suggested by Ellison et al. (|2008h . 
who presented a mild (~ 0.1 dex) dependence of metallicity 
on SSFR when binning galaxies a ccordi ng to their stellar 
mass. After MIO, Lara-Lopez et al. l|2010l ) discussed a plane 
in the same 3D space, and proposed to use it to derive the 
stellar mass of galaxies once their SFR and metallicity is 
known. The exact shape of the relation, especially at the 
high mass end, depends on the choice of the metallicity 
calibrations, on the assumptions about dust extinction and 
SFR, and on the extent of the ap erture corrections allowed, 
as reported in Yates et al. (|201ll ). 

Although the presence of the FMR is well established 
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at low z, regardless of the metallicity calibration used, only 
a limited amount of distant galaxies have been used for 
the study of the FMR at high redshifts. The data available 
show that no evolution is present up to z=2.5, see Fig. 4 of 
MIO), i.e., galaxies of any mass and SFR follow the same 
FMR defined by local galaxies. The observed evolu tion o f 
the mass-metallic ity re lation (e.g. Maiolino et al. (|2008h . 
Mannucci et al. (|2009l) ) is therefore due to a prospective 
effect, as high redshift galaxies usually have higher SFR 
and, as a consequence, lower metallicities. Some evolution 
is observed at even higher redshifts, z>3, as discussed in 
MIO and Troncoso et al. (in preparation). 
This results has been later c onfirm ed by several groups. 
For example, Richard et al. l|201ll ) studied a sample of 
lensed galaxies at 1.4<z<5, which have SFR smaller than 
the other s ample s of un-lensed gal axies at these redshifts 
(Erb et al |2006| . MaioUno et al. |2008| . Mannucci et al. 
I2OO9I ). These galaxies show higher metallicity with respect 
to the mass-metallicity relation, exa ctly a s expected by 
the FMR. Similarly, Nakajima et al. (|201ll ') presented the 
metallicities of stacked Lya emitter at z ^ 2.2, which have 
smaller SFR than other galaxies at similar redshift, but 
also higher metallicities according to the FMR. On the 
other hand, various highly star forming galaxies at high 
redshift are presented in the literature as very metal poor 
for their stellar mass. However, they result to be consistent 
with t he FM R once their h igh S FR is considered (e.g. Erb 
et al. l20ld . Kassin et al. boili). Fi nallv. the use of the 
FMR allowed Mannucci et al. ( 201ll ') to show that, when 
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considering their high SFRs, the host galaxies of long-GRBs 
have the same metallicity distribution of field galaxies. 

The high redshift galaxy samples used by MIO show a 
much larger spread in metallicity (~0.25 dex) around the 
FMR than local SDSS galaxies (~0.05 dex). The spread is 
comparable with the much larger uncertainties in SFR and 
metallicity, due to the generally lower signal-to-noise (S/N) 
ratio of the detected emission lines. This large spread could 
hide systematic trends, therefore it is important to test this 
results using larger samples of galaxies. In this paper we use 
data from the zCOSMOS survey to study the evolution of 
the FMR up to z ~ 0.8 with higher precision. This is funda- 
mental to understand if the same long lasting equilibrium 
between star formation, infall of pristine gas and outflow 
of enriched material that seems to be res ponsib le for the 
F MR in local galaxies (see e.g. Dave et al. I2OIII . Torrey et 
al. l201lh is also driving galaxy evolution at higher redshift, 
wher e infaUs (Erb I2008i . Mannucci et aL 12009 . Cresci et al. 
I2OI0I) and outflows (e.g. Steidel et al. |201(F are expected to 
be even more relevant. 

The paper is organized as follows. In Sect. [2] we briefly 
introduce the sample used, in Sect. |3] we present the method 
used to derive the stellar masses and in Sect. S] we described 
how the line fluxes and the flnal sample selection were ob- 
tained. The measures of SFRs and extinctions are discussed 
in Sect. [SI while those of metallicities in Sect.[51 Our results 
are discussed in Sect. [8] while our conclusions follow in Sect. 

m 



2 THE ZCOSMOS SAMPLE 

zCOSMOS (Lilly et al. |2007| ) is a redshift survey in the 
COSMOS field using the VIMOS multi-object spectro- 
graph on ESO/VLT. It consists of two parts: the first, 
zCOSMOS-bright, is a redshift survey of - 20000 I-band 
selected galaxies at redshifts z < 1.2, covering the approx- 
imately 1.7 deg^ of the COSMOS field; the second part, 
zCOSMOS-deep, is observing ~ 10000 galaxies selected 
through well-defined color selection criteria and at Bj25 to 
be in the redshift range 1.5 < z < 3.0. These observations 
are restricted to the central 1 deg^ of the COSMOS field. 
Here we have used the spectra of the data release 2 of 
the zCOSMOS-bright survey, the so-called "lOk-bright" , 
available on the project web site since 1 Oct 200fl These 
consist of ~ 10000 spectra, obtained with 1" wide slits and 
the MR grism (R~600 with a wavelength range of approx- 
imately 5550 to 9650 A, sampled at roughly 2.5 A/pixel). 
The details about the object selection, o bserv ations and 
data reduction can be found in Lilly et al. (|2009l ). 
In order to measure gas phase metallicity for our galax- 
ies using strong optical emission lines, we need either 
[NII]A6584 and Ha, or [OII]A3727, [OIII]A5007 and H/3 
visible in the wavelength range of the VIMOS spectrum. We 
therefore divided the obtained sample in two subsample at 
different redshifts: 2825 galaxies at z < 0.47 have accessible 

^ http:/ /archive. eso.org/cms/eso-data/data-packages/zcosmos- 
data-release-dr2/ 



Ha emission and 3691 galaxies at 2 > 0.49 are in the right 
redshift range to observe [OH] emission. [OIII]A5007and 
H/? are observable for all these galaxies. 



3 MASS AND CONTINUUM SPECTRA FROM 
SED FITTING 

To estimate the stellar masses we used 12 photometric 
bands: CFHT u* and Ks, Subaru Bj, Vj, g+ , r+ , i+ and 
z+, UKIRT J (Capak et al. 2007, Ilbert et al. 2009), and 
Spitz er IRAC at 3.6 ^m, 4.5 ^m and 5.8 ^m (Frayer et al. 
[2009h . 

The photometry for all the optical and near-IR bands were 
measured over an aperture of 3" diameter, to minimize 
the effect of Point Spread Function variation from band to 
band. The derived fluxes are then corrected to total using 
the so called "auto-offsets", derived from the i-band image 
comparing the 3" magn itude with the SExtractor total one 
(see Capak et al. \200t ). The IRAC fluxes were measured 
in a circular aperture of radius 1.9", converted to total 
fluxes using the correction factors derived by Surace et al. 
(|2004[ ). We also used G AL EX FUV and NUV photometry 
from Zamojski et al. |2003), when available. A few galaxies 
show problematic values of the magnitude (Bolzonella et 
al. private communication), therefore GALEX photometry 
was only used when the reduced of SED fitting is below 
2, for a total of 64 galaxies. 

The stellar masses were obtained with a Spectral 
Energy Distribution (SED) fitting ana lysis, using the 
Hype rZmass code (Pozzetti et al. |2007| . Pozzetti et al. 
l2009l ). a m odified version of the HiperZ code (Bolzonella 
et al. I2OO0I I : the observed photometric SEDs are compared 
to those obtained from a set of reference synthetic spectra 
from stellar population models, and the best fitting SED 
is obtained using a minimization. The best fitting 
spectrum and its normalization provide estimates of the 
star formation rate, extinction, stellar population age and 
galaxy stellar mass (as the integral of the SFR over the 
obtained age, correcting for mas s losses during stellar 
evolution, see Pozzetti et al. l2007h . In addition, the code 
makes available in output the best fitting synthetic stellar 
continuum for each galaxy, that has been used to correct 
for the Balmer absorptions (see Section 3)). 
The refere nce sp ectral library was built with the Bruzual 
& Chariot l|2003l ) code for spectral synthesis models, using 
the low resolution version with the Padova 1994 tracks. 
The models were produced assuming a Chabrier Initial 
Mass Function (IMF) (Chabrier et al |2003i ) with an upper 
mass limit of IOOM0, the Calzetti |200(7) extinction law 
with Av varying between — 4 mag, and fixing the redshift 
to the spectroscopic value. We used smooth exponentially 
decreasing Star Formation Histories, with time scale r in 
the range 0.1 — cxd and age t between 0.1 — 20 Gyr. An 
example of SED fitting for one of the galaxies of our sample 
is shown in Fig. [1] 

The stellar masses computed by Kauffmann et al. (|2003ah 
for the SDSS sample used by MIO to derive the FMR in 
the local universe were instead computed using random 
secondary bursts in addition to the single one used here. 
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Figure 1. Example of SED fitting for one of the galaxies in the 
sample, ZC803435. The observed photometric data are shown as 
black solid points, while the expectations for each band based on 
the synthetic spectrum are shown as empty points. 



Pr evious studies (e.g. Fontana et al. |2004| . Pozzetti et 
al. l2007h have shown that random secondary bursts may 
produce systematically larger masses up to about 0.2 dex. 
We therefore run several tests adding random secondary 
burst to our synthetic spectra to check if this systematic 
difference in the stellar mass due to the use of a single 
exponential star formation history in our data is present 
and comparable to previous studies. In agreement to what 
is found in the literature, we found that the stellar mass 
obtained with random secondary burst are ~ 0.15 — 0.2 
dex larger than the ones obtained with a single burst. We 
therefore added 0.2 dex to the stellar masses obtained with 
a single burst for consistency with the SDSS stellar masses. 
The stellar masses derived were also com pared with the 
results obtained using the Maraston et al. (|2005l ) spectral 
library, which include a more detailed treatment of the 
thermally pulsing asymptotic giant branch stars. These 
stars contribute significantly to the total luminosity in the 
~ 0.2 — 2 Gyr stellar population age. We found th at the 
stellar masses obtained using the Maraston et al. (|2005l ) 
library are lower on average by a factor of ~ 0.1, as already 
mentioned in Bolzonella et al. \2Q1(\ ). To allow an easy 
comparison with the previous study in MIO for SDSS 
galaxies, we use in the follow ing the results obtained with 
the Bruzual & Chariot (|2003l ') library. 



4 LINE FLUXES 

The extracted 1-D spectra provided by the zCOSMOS 
release were used to measure emission line fluxes. For 
each spectrum we fitted several emission lines, [OIIJA3727, 
[NeIII]A3869, H7, H/3, [OIII]A5007, [OIII]A4959, [NII]A6548, 
Hq, [NII]A6584, [SII]A6716, [SII]A6730, when available in 




Figure 2. Example of the measurement of line fluxes by using 
the continuum obtained with SED fitting for one of the galaxies 
in the sample, ZC817676. The observed spectra is shown in black, 
while the continuum model and the best fitting Gaussian for the 
H/3 line are shown in red. The [OIII] doublet is also visible in the 
spectral region shown. 



the accessible spectral range (see Section [Sjl . The fit was 
performed using a custom IDL code, which automatically 
fits with a minimization a suitable stellar continuum 
and a Gaussian to each emission line. 

To measure reliable fluxes, it is important to properly 
account for the stellar Balmer absorption, which in some 
galaxies can reach equivalent widths of ^ SA. The best 
fitting synthetic stellar continuum derived from the SED 
fitting analysis of each galaxy (see Section [Jl is used to 
reproduce the absorption line profile and the continuum 
shape. The synthetic stellar continuum is providing a 
very good representation of the observed continuum in 
all the inspected spectra, providing a good estimate of 
the absorption underlying the Balmer emission lines: an 
example is shown in Fig. [51 

The continuum is fitted in a region of ±300 A around 
the expected position of each emission line. During the 
minimization, pixels in the data that consistently deviated 
more than 3cr from the average were rejected, and were not 
used in the analysis. The uncertainties on the measured 
centroid position, velocity dispersion and fiux are evaluated 
through Monte Carlo realizations, perturbing the input 
data assuming Gaussian uncertainties. 



5 STAR FORMATION RATE AND 
EXTINCTION 

SFR can be derived from the Ha and H/3 lines, which pro- 
vide a direct probe of the young massive stellar population. 
Similar to what has been done in MIO and many other 
authors, w e hav e used the classical conversion factor by 
Kennicutt l|l998t ). scalin g down the results by a factor of 
1 /^(Pozzetti et al. |2007| ) to convert them to the Chabrier 
((20031 ) IMF. The conversion factor is based on the Ha or 
H/3 flux corrected for dust extinction, and therefore requires 
a good knowledge of the amount of extinction Av suffered 
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by the emission lines. 

For galaxies at z < 0.47, where both Ha and H/3 are visible, 
the extinction can be derived directly from the Balmer 
decrement Ha/H/3, assuming a Calzetti (,2000. ) extinction 
law, and the SFR from the Ha extinction corrected flux. At 
z > 0.49, where Ha is not visible anymore in our data, we 
used the Balmer decrement of the fainter H7/H/3 for all the 
galaxies where S/N(H7)J^ 5, while the extinction Av,sed 
derived for the continuum by the SED fitting (see Section 
[si for the other galax ies. 

Calzetti et al. (I2OO0I ) found evidence that the extinction 
derived from emission lines Av,neb in actively star forming 
galaxies is not equal but proportional to that measured 
from the continuum, with Av,neb_ = Av,sbd/0.44 (see also 
e.g. Forster Schreiber et al. |2009| ). We have compared 
Av,neb and Av,SED for all the galaxies where a Balmer 
decrement extinction measurement is available (either 
Ha/H/3 or H7/H/3), finding that the extra attenuation of 
the emission lines is confirmed by our data: our best fitting 
relation is Av.se d = (0.42 ± 0.01) ■ Av.net- Given the small 
differences, we adopt the relation Av,neb = Av,sed /O-'l'i of 
Calzetti et al. (200(1 ) to all the objects in the higher redshift 
bin with no H7 flux, in order to ensure a direct comparison 
of our results with previous works. 



6 METALLICITY 

Given the large variety of different strong line metallicity 
calib rations used in the l iterature (see e.g. Pagel et al. 
I1979I . Pettini & Pagel |2004| : Ke wley fc Dopita l2002l : Dopit a 
et al. I2OO6I : Tremonti et al. |2004| : Nagao et al. |2006| ). 
which give very different results for a given line ratio, it 
is mandatory to study the evolution of the FMR using 
the same calibration used in MIO. Consistently with that 
work, gas-phase metallicities were c ompu ted using the 
calibr ations obtained by Nagao et al. (|2006l ) and Maiolino 
et al. (|2008l l. to allow a proper comparison with local SDSS 
galaxies. These ratios have been empirically calibrated at 
low metallicities (12 + log{0/H) < 8.4) using direct Te 
measurements in local metal poor galaxies and HII regions, 
and are based on the photoionization models by Kewley 
& Dopita lioO^ (hereafter KD02) at high metallicity. The 
metallicity is therefore provided on the KD02 scale, while 
a conversion to other met allicity systems can be found 
in Kewley & Ellison (|2008l '). The calibration proposed in 
KD02 is dependent on the [Nil] /[OH] ratio, but the use 
of a secondary element as Nitrogen is of some concern. 
This issue will be addressed in detail by Maiolino et al . 
(in preparation, see also P erez-Montero & Contini l2009al . 
Torres-Papaqui et al. I2OIII ). 

Two different measurements of metallicity are available 
for the two redshift ranges. At 2: < 0.47, where Ha is still 
visible, the [NII]/Ha ratio is used. At z > 0.49, where 
[OH] emission is in the accessible wavelength range, the 
measurement is instead based on the R23 quantity, deflned 
as 7l23=([OH]A3727+[OIII]A4959+[OHI]A5007)/H^. In this 
latter the lines involved are significantly separated 

in wavelength, the extinction correction discussed is Section 
[5] has to be applied to the observed line fiuxes. Moreover, 



as R23 has two possible metallicity solutions for a given 
line ratio, other diagnostics have to be used to remove the 
degeneracy and distinguish between the two branches at 
low and high metallicity. 

We use for this task both the [OIII]/[OII] and 
[NeIII]/[OII] line ratios, which have a monotonic depen- 
dence on metallicity, but are also affected by a larger 
dispersion due to their additional d epend ence on the 
ionization parameter (see Nagao et al. I2OO6I ). In spite of 
their larger scatter, low and high metallicity galaxies are 
reasonably well s eparat ed in both diagrams (see e.g. Fig. 5 
in Maiolino et al. l2008l '). 

We therefore assigned all galaxies with log([OIII]/[OII])< 
0.45 and log([NeIII]/[OII])< -0.6 to the high metallicity 
R23 branch, while all the others to the low metallicity 
one. The two independent branch indicators used agree for 
all the galaxies in our sample, supporting our method to 
discriminate between branches. 

To further check the reliability of the branch sele ction done, 
we add to the galaxies used in Nagao et al. (|2006l ) and 
Maiolino et al. (2008) the sample of ~ 6 00 SDSS with direct 
Te measurement of Kniazev et al. (|2004h . We find that their 
galaxies with S/N > 5 on the [OII I]A43 63 line (used for Te 
measurement, see e.g. Izotov et al. |2006| '1 follow our criteria 
on [OIII]/[OII] to discriminate between the two branches. 
This is also confirmed by a larger sample of galaxies with 
Te measurement that will be presented in a forthcoming 
paper (Maiolino et al. in preparation). 

The uncertainties on the gas phase metallicities are of 
the order of 0.1 — 0.2 dex, and are dominated by the spread 
in the calibrations themselves. The uncertainties in the 
higher redshift bin tend to be larger when the metallicity 
has values 12 -I- log(0/H) ~ 8, where R23 has a maximum 
and a small uncertainty in the line ratio produces a large 
uncertainty in metallicity. In any case, we note that just 
very few galaxies in our sample are in the low metallicity 
region where R23 is fiat (between 12-|-log(0/H)'-^ 7.7 — 8.4, 
see Fig. g)) 



7 SAMPLE SELECTION 

The results of the SED and emission line fitting described 
in Sect.[3]and[4]above were combined and used to define the 
two final subsamples (galaxies at z < 0.47 with accessible 
Ha emission and galaxies at, z > 0.49 with [OH] emission). 
We applied to the initial sample discussed in Section [2] the 
following selections: 

(i) each zCOSMOS spectrum is assigned to a confidence 
class according on the reliability of the final redshift. We 
selected all the objects with confidence classes 3 and 4 (re- 
spectively, "secure" and "very secure" redshifts), and class 
2.5 (redshift with reliability larger than 85% and photomet- 
ric redshift in good agreement with the spectroscopic red- 
shift). We also added all the objects in classes 23 and 24, 
i.e. secondary targets detected in a slit positioned on an- 
other target, but with secure and very secure redshifts. We 
further excluded from our sample duplicate objects. Galax- 
ies classified as broad-line AGN in the zCOSMOS catalog 



The metallicity properties of zCOSMOS galaxies at 0.2<z<0.8 5 



(class > 10) or stars (class *.4) are automatically excluded 
from our selection. 

(ii) as our two subsamples may be still contaminated by 
type-2 narrow line AGNs, or with by AGN with nuclear 
emission diluted in the host galaxy light, we applied two 
additional filterings. First, we excluded objects associated 
to X-rays selected AGN from the XMM-COSMOS cata- 
log with rest-frame hard X-rays luminosity in ex cess th an 
2 • 10*^ erg s"\ typical of AGNs (see Brusa et al. I2OI0I for 
more details). Second, to further exclude Compton Thick 
and Low Luminosity AGN below the XMM-COSMOS sen- 
sitivity, we used standard BPT selection s base d on emission 
line ratios. Following Bongiorno et al. (|2010l ). we used the 

[oiir 

et al. 



<z> = 0.25 



<z> = 0.60 



/H/3 v ersus [Nil] /Ha diagnostic diagram (Kauffmann 
2003bh for galaxies &t z < 0.46, and the [O III]/H/3 ver- 



sus [OII]/H/3 diagram (Lamareille et al. l2004f ) for galaxies 
at 2 > 0.49; 

(iii) the velocity shift between each fitted emission line 
centroid and the zCOSMOS redshift must be below 280 
km/sec, and the line width must be between 110 km/sec 
and 290 km/sec, otherwise the fit is considered unreliable 
and the galaxy is excluded. For this selection we consid- 
ered all the main emission lines: H/? and [OIII]A5007 for all 
the galaxies. Ha and [Nn]A6548 for galaxies at low redshift 
(z<0.46), [OII]A3727 for galaxies at high redshift (z>0.49); 

(iv) the same lines must not correspond to bright sky 
lines, or to spectral regions badly affected by fringing (A > 
8000 A); 

(v) the reduced of SED fitting must be below 5, and 
the final 68% uncertainty on mass must be below 0.30 dex; 

(vi) a threshold on the S/N of Ha and H/3 is also applied 
to have reliable and unbiased metallicity measurements. As 
line ratios depend stron gly on metallicity (see, for exam- 
ple, Maiolino et al. l2008l l. applying the same thresholds to 
all relevant lines or to their ratios would introduce severe 
biases in metallicity (see MIO for a detailed discussion). 
To avoid these biases, and following what has been done 
for MIO, we only apply high S/N thresholds to the lines 
which are diagnostics of SFR like Ha and H/3. The actual 
threshold values depend on the redshift and stellar mass: 
at z<0.46 we used S/N(Ha) = 15 for log(M^/MQ) ^ 9.4, 
and S/N(Ha)=30 for log{M»/MQ) < 9.4. At z > 0.49 we 
used S/N(H/3)=8 for log{M:>/MQ) 9.8, and S/N(H/?)=20 
for log{Mt/MQ) > 9.8. With these choices aU the other 
relevant lines are detected with adequate S/N without ad- 
ditional, metal dependent selections. Indeed, in our lower 
redshift bin biases on metallicity are important mainly at 
low stellar masses: for example, at log{M / Mq) ~ 9, where 
metallicities are 12-|-log(0/H)-- 8.5, [Nil] is only ~ 0.1 Ha, 
while at log{M / Mq) ~ 10.2, where the metallicity is as high 
as 12-|-log(0/H)~ 9, [Nil] is ~ 0.4 Ha, and a lower S/N 
threshold on Ha may be applied to maximize the number 
of galaxies in our sample. In the higher redshift bin, on the 
other hand, the situation is reversed and [OIII], which is 
only ~ 0.4 H/3 in the high mass bin, increases to ~ 3 H/? at 
the lower masses and metallicities. 

These considerations indicate that our choice of S /N thresh- 
olds on Ha and H/3 lines that depends on mass in both red- 
shift bins maximize the number of galaxies selected without 
introducing metallicity biases; 

(vii) following MIO, we checked the metallicity derived 
from [NII]A6548/Ha for galaxies at 2 < 0.46 using an in- 
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Figure 3. Redshift distribution of the two final galaxy subsam- 
ples (see text). The galaxies at 2 < 0.47 with visible Ha emission 
are shown in red, and galaxies at z > 0.49 with [Oil] emission are 
shown in blue. 



dependent metallicity indicator. As [OH] is not available in 
this redshift range, we used [0III]/H/3, which has proper- 
ties similar to R23 with just a slightly larger scatter: we se- 
lected only galaxies where the values of metallicity obtained 
with the two different indicators differ less that 0.25 dex. 
At 2 > 0.49, instead, no additional independent metallicity 
tracer with low enough scatter is available, and no additional 
selection is applied; 

After applying all these selections, we obtain 169 galax;- 
ies at 2 < 0.47 (< 2 >= 0.25) and 165 galaxies at 2 > 0.49 
(< 2 >= 0.60). The redshift distribution of the two final 
subsample is shown in Fig. [3] 



8 RESULTS 

The Mass-Metallicity relations for the two zCOSMOS 
samples at 2 ~ 0.25 and 2 ~ 0.60 are shown in Fig. U (left 
panels). As expected, the Mass-Metallicity relation for the 
z ~ 0.25 sample is comparable to the relation derived by 
MIO for SDSS galaxies, as the SDSS sample used by MIO 
is including galaxies up to 2 < 0.3. 

In the higher redshift bin, consistently with what was 
observed in other samples at compa rable redshift (e.g. 
Savaglio et al. I2005i . Co wie fc Barger I2OO8I . Rodri gues et 
al. 200i, La marei lle et al. |2009| . Perez- Motero et al. ,2009b|, 
Zahid et al. l201ll ). we find significant evolution of the mass 
metallicity relation, of ~ 0.1 dex. Moreover, if we plot the 
single galaxies on the mass metallicity relation color coded 
as a function of their SFR, as shown in Fig. [S] (left panels) 
for both redshift bins, there is a trend for more star forming 
galaxies to have lower metallicity at a given stellar mass, as 
found by MIO for the SDSS sample. 

MIO, in fact, defined the FMR adding the SFR of the 
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Figure 4. The mass metallicity relation (left panel) and the /io.32 projection of the Fundamental Metallicity Relation (right panel) for 
the 2 < 0.46 (red solid lines) and the z > 0.49 (blue solid lines) zCOSMOS galaxies. The shaded areas contain galaxies between the 16th 
and the 84th percentile rank. The two high redshift relations are compared with the local one, as defined by the SDSS galaxies from 
MIO (black solid lines). It is evident an evolution of the mass metallicity relation at 2 > 0.49, while the corresponding projections on the 
FMR plane shows no evidence of evolution even at the higher redshift. 



galajdes as a third parameter in the mass metallicity rela- 
tion, showing that SDSS galaxies with lower metallicity also 
have higher SFR at a given mass, and that roughly half of 
the scatter in the relation was due to the different individual 
SFR of the galaxies. They also showed that the metallicity 
scatter of local galaxies could be minimized plotting them 
on a particular projection of the FMR, introducing a new 
quantity /io.32 defined as: 

= log{M/MQ) ~ 0.32 ■ log{SFR) (1) 

As shown in the right panel of Fig. [H no evolution is 
observed even in the higher redshift bin in the /io.32 plane, 
once the increased SFR of the zCOSMOS galaxies at high 
z with respect to the local SDSS galaxies is accounted for. 
Moreover, the projection clearly removes any dependence 
of the mass-metallicity scatter from the SFR of the galaxies 
(see Fig. O right panels). This shows that the star forming 
galaxies up to z ~ 0.8 populate the same tight surface in 
the 3D space defined by their mass, metallicity and SFR as 
demonstrated by MIO for local SDSS galaxies. 

To better check the presence of any evolution of the 
FMR with redshift, we compute the expected metallicity of 
each galaxy, given its stellar mass and SFR, according to 
the 3D plane defined by local SDSS in MIO. We then plot in 
Fig. [6] the difference of the observed value with the expected 



metallicity. The galaxies in both redshift bins are consistent 
with no differences from the FMR, with an average offset of 
9 ■ 10"^ ± 0.17 at z < 0.46 and 9 • 10"=^ ± 0.16 at z > 0.49, 
confirming with a much larger sample the no evolution of 
the FMR up to z ~ 0.8 suggested by MIO. 

The larger observed scatter in the zCOSMOS data 
with respect to the SDSS galaxies of MIO is partly due 
to instrumental effects. In fact, the VIMOS spectra used 
here are affected by significant fringing beyond ~ 8000 
A (see e.g. Le Fevre et al. 120051 ). Although we excluded 
galaxies with emission lines in particularly bad part of 
the spectrum, our line fitting quality is anyway degraded 
redward of ~ 8000 A. In fact, if we exclude the galaxies 
with emission lines contaminated by fringing (i.e. galaxies 
at z > 0.22 for our lower redshift bin and objects at z > 0.6 
for the higher redshift bin), the scatter is significantly 
reduced to 0.11 dex and 0.14 dex respectively for the two 
subsamples. Moreover, the scatter in the local FMR based 
on SDSS galaxies depends significantly on mass and SFR, 
increasing at higher SFR and lower stellar mass (see Table 
1 in MIO). Therefore, the scatter derived for the zCOSMOS 
objects has to be compared with that obtained for SDSS 
galaxies with comparable stellar mass and SFR only: the 
z — scatter corresponding to the lower redshift zCOSMOS 
subsample galaxy properties is 0.10 dex, while the one 
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Figure 5. The mass metallicity relation (left panels) and the 
/iO.32 projection of the FMR (right panels) for the individual 
galaxies of the sample. The lower redshift bin {z < 0.46) is shown 
in the upper panels, while the higher redshift bin (z > 0.49) in 
the lower panels. The galaxies are color coded as a function of 
their SFR, with decreasing SFR from blue to magenta. While a 
trend for lower metallicity at higher SFR and fixed mass is present 
in the mass metallicity relation at both redshifts, the projection 
removes such effects, showing that the galaxies up to z ~ 0.8 
populate a 3D plane defined by mass, metallicity and SFR as 
local SDSS galaxies. 

relative to higher redshift one is 0.11 dex, closer to the 
fringing free scatter observed in our data. The remaining 
difference in the higher redshift bin is probably due to the 
single metallicity indicator used for those galaxies, while in 
MIO the abundances were derived using both R23 and Hq 
for all the galaxies. 

The redshift evolution of the FMR based on the zCOS- 
MOS data presented in this work and on higher redshift data 
from the literature compiled by MIO is shown in Fig. [T] All 
the galaxy samples up to z=2.5 are consistent with no evo- 
lution of the FMR defined locally, while metallicities lower 
by ~ 0.6 dex are observed a.t z — 3.3 (see MIO). As the data 
points at 2: > 1 are based on few massive galaxies, the larger 
samples at high redshift presented here at 2: ~ 0.25 and 
z ~ 0.63 provide a stronger evidence for the no-evolution of 
the FMR with redshift. 



9 CONCLUSIONS 

In this paper we have studied the gas phase metallicity of 
a sample of more than 300 star forming galaxies between 
0-2 < z < 0.8 from the zCOSMOS survey. We have mea- 
sured the gas phase metallicity using strong optical line ra- 
tios, [NII]A6548/Ha for galaxies at 2: < 0.46 and i?23 for 
galaxies at z > 0.49. 

We find that the metallicity of the galaxies is tightly re- 
lated to their stellar mass, as for local galaxies. However, 



Figure 6. Difference from the FMR surface for the zCOSMOS 
galaxies at z < 0.46 (left panel) and 2 > 0.49 (right panel) as a 
function of stellar mass. The single points are shown color coded 
as a function of their specific SFR, while their average in each 
bin along with the la scatter is shown in black. The galaxies 
in both redshift bins are consistent with no differences from the 
FMR, with an average ofl'set of 9 • 10~^ ± 0.17 at 2 < 0.46 and 
9 • 10~^ ± 0.16 at 2 > 0.49. No residual dependencies from the 
SSFR is found. 



a clear evolution of the mass metallicity relation of ~ 0.1 
dex with respect to the local one is observed for the galaxies 
with z > 0.49. We also found a dependence of the scatter 
observed around the mass metallicity of our sample on the 
SFR, showing that star forming galaxies up to 2 ~ 0.8 are 
defining a 3D plane between stellar mass, metallicity and 
SFR, as local SDSS galaxies. Moreover, we have shown that 
projecting the galaxies on the FMR plane using the po.32 
quantity successfully removes the dependence of the scatter 
by the SFR of the galaxies, and that on this projection the 
higher redshift zCOSMOS galaxies follow the relation de- 
fined by local SDSS galaxies with no evidences of evolution. 
The lack of evolution of the FMR at 0.2 < z < 0.8 is fully 
confirmed for the zCOSMOS sample, as the expected metal- 
licities computed according to the local FMR using the mass 
and SFR of the galaxies are consistent with the observed 
ones: the average offset is 9 • 10~^ ± 0.17 at z < 0.46 and 
9 ■ 10~^ ± 0.16 at z > 0.49. The observed scatter is somehow 
larger than the one measured for SDSS galaxies, although 
the difference is probably due to VIMOS instrumental ef- 
fects and to the lower number of independent metallicity 
indicators available. 

The results presented confirm with much higher level of sig- 
nificance that the same physical processes are probably re- 
sponsible for the SFR, metal enrichment and gas recycling 
through infall and outflows in the local and higher red- 
shift Universe, as already suggested by MIO. Moreover, they 
support the idea that the observed evolution of the mass- 
metallicity relation, at least up to z ~ 0.8, is simply due to 
more severe selection biases at high redshift as well as to 
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Figure 7. Redshift evolution from the FMR plane defined by 
local SDSS galaxies for zCOSMOS galaxies (this work, red and 
blue) and other high-z sample in the literature (0.5 < z < 0.9: 
Savaglio et al. 2005; 1.0 < 2 < 1.6: Shapley et al. 2005, Liu et al. 
2008, Wright et al. 2009, Epinat et al. 2009; 2.0 < z < 2.5: Law 
et al. 2009, Forster Schreiber et al. 2009; z ~ 3.3: Maiolino et al. 
2008, Mannucci et al. 2009). All the galaxy samples up to z=2.5 
are consistent with no evolution of the FMR defined locally, while 
metallicities lower by ~ 0.6 dex are observed at 2 = 3.3. While 
each black point at 2 > 0.7 is the average of ~ 20 galaxies or less, 
the points at 2 ~ 0.25 and z ~ 0.60 comprise ~ 160 galaxies each. 



the increase of the average SFR, which results in samphng 
different parts of the same FMR at different cosmic epochs. 
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